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ABSTRACT 

We propose a test that can in principle detect any systematic errors in the Hipparcos 
parallaxes towards the Hyades cluster at the level of 0.3 mas. We show that the statistical 
parallax algorithm subsumes the classical moving cluster methods and provides more 
precise estimates of the distance and the hrst two moments of the velocity distribution of 
the Hyades cluster namely, its bulk space velocity and the velocity dispersion tensor. To 
test the Hipparcos parallaxes, we hrst rescale the bulk velocity determined from statistical 
parallax to force agreement with the distance scale determined from Hipparcos parallaxes. 
We then predict the parallaxes of Hyades cluster members using this common cluster space 
velocity and their Hipparcos proper motions. We show that the parallaxes determined in 
this manner (vTpm) are consistent at the la level with the parallaxes (TTorb) of three Hyades 
spectroscopic binary systems with orbital solutions. We hnd that (TTpm — vTorb) = 0.52 ±0.47 
mas, where the error is dominated by the errors in the orbital parallaxes. A reduction in 
these errors would allow a test of the systematic errors in the Hipparcos parallaxes at the 
0.3 mas level. If the Hyades distance scale is hxed by Hipparcos parallaxes, then its bulk 
velocity in equatorial coordinates is (14, Vy, Y) = (—5.70 ± 0.20,45.62 ± 0.11, 5.65 ± 0.08) 
kms“\ its velocity dispersion is 320 ± 39 ms“^, and the distance modulus to the centroid 
of our sample of 43 cluster members is 3.34 ± 0.02 mag. 

Subject headings: astrometry: parallaxes, methods: analytical, statistical. 

Galaxy: open clusters and associations: Individual(Hyades) 


1. INTRODUCTION 


The Hipparcos mission ([ESA97|) has derived accurate astrometric parameters for about 
120,000 stars distributed all over the sky. The systematic errors in the Hipparcos parallaxes 
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are estimated to be ^ 0.1 mas, while the random errors are of the order of 1 mas 


renou 


et al. 19^ [Arenou, Mignard fc Palasi 1997| ). However, recent comparisons of the distances 
to open clusters derived from Hipparcos parallaxes and main sequence htting techniques 
show surprisingly large differences for some clusters ( IVlermilliod et ah 1997| ; Ho bichon et 
ah 1997[) , which can be reconciled if the systematic error in the Hipparcos parallaxes is 
at the level of 1 mas, at least on small angular scales (Pinsonneault et al. 1998, hereafter 
PSSKH98). It is therefore prudent to compare the Hipparcos parallaxes with accurate 
parallaxes determined in an independent manner. In this paper, we propose and apply one 
such method which could in principle test for systematic errors in Hipparcos astrometry 
towards the Hyades cluster at the level of 0.28 mas. 


The statistical power of our test arises from the fact that the fractional error in the 
Hipparcos proper motions of the Hyades cluster members (a^/ (/iHya) = 1.4%) is about four 
times smaller than the fractional error in their Hipparcos parallaxes (cijr/ (vrHya) = 6%), 
where we have assumed that the mean parallax and the proper motion of the Hyades cluster 
are (vTHya) = 21.5 mas and (/^Hya) = HI masyr“^ respectively and their errors are cx^ = 1.3 
mas and = 1.5 rnasyr”^. Hence, if one can accurately determine the space velocity of the 
Hyades cluster, one can use the Hipparcos proper motions to predict the parallaxes of the 
individual Hyades members more accurately than the Hipparcos parallaxes, assuming that 
all the members partake in this common cluster motion (to within the velocity dispersion 
of the cluster). The accurate parallaxes predicted in this manner can then be compared 
with the Hipparcos parallaxes and with parallaxes determined in an independent manner 
to determine the level of the systematic errors in the Hipparcos astrometry. 


The basic steps of our test are as follows: 


(1) : We select a set of Hyades cluster members all of which are consistent with having the 

same velocity as the bulk motion of the cluster itself. 

(2) : We derive the space velocity of the Hyades cluster by combining the radial velocities of 

the individual cluster members, their Hipparcos proper motions and their photometric 
distances (up to an initially unknown global scale factor) derived from an isochrone of 
the Hyades main sequence. All these independent data are elegantly combined in the 
statistical parallax method to derive a maximum likelihood estimate of the distance 
scale and the hrst two moments of the velocity distribution of the Hyades - the bulk 
velocity and the velocity dispersion tensor ([Hawley et al. 19"^ ; jStrugnell, Reid fc 
Murray 1986| ; Popowski fc Gould 1998]) . This procedure generalizes the classical 
moving cluster methods. 


(3): We rescale the distance scale parameter of the statistical parallax solution from (2) 
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to force agreement with the Hipparcos trigonometric parallaxes of the member stars 
in (1). This yields the best estimate of the bulk velocity of the cluster under the 
assumption that Hippareos parallaxes do not suffer from systematic errors. 

(4) : We adopt this common cluster velocity to predict the parallaxes of three Hyades 

spectroscopic binary systems whose orbital solutions are known accurately from 
previous work. 

(5) : We compare the parallaxes predicted in this manner with the parallaxes from the 

orbital solutions of the three binaries to check if there are any systematic errors in the 
Hipparcos astrometry towards the Hyades cluster. 

The outline of this paper is as follows. We describe the various assumptions underlying 
our method in §2. We explain the connection between the classical moving cluster methods 
and the statistical parallax method in §3. In §4, we describe our selection of Hyades cluster 
members and use the statistical parallax algorithm to determine the common space velocity 
and the distance to the centroid of the cluster. We outline the procedure for estimating 
the parallaxes to individual cluster members using the common cluster motion and proper 
motions from the Hipparcos catalog in §5. We compare the parallaxes predicted by this 
method with the parallaxes from the orbital solutions of the 3 binary systems in §6. We 
summarize our conclusions in §7 and describe the future potential of this technique. As an 
aside, we note that we will drop the usual conversion factor Ay = 4.74047 kmyrs“^ from all 
our equations for the sake of clarity, leaving it to the reader to include it in the appropriate 
equations. 


2. ASSUMPTIONS 

We now describe the assumptions underlying the technique described in this paper. 
These include assumptions about both the kinematical structure of the cluster and the 
quality of the observational data. These assumptions are: 

(1): The velocity distribution of the Hyades cluster members can be adequately described 
by an average bulk velocity and an isotropic intrinsic velocity dispersion tensor. In 
particular, we assume that the velocity structure of the cluster is not signihcantly 
affected by shear, rotation or expansion. From the Hyades age of T 625 Myr ( TO) , 
its velocity dispersion, a ~ 330 ms“^, and its characteristic size a ~ 5 pc, we know 
that the cluster has survived for aT /a ~ 40 crossing times. Hence, it is gravitationally 
bound and it is likely that all non-conserved modes of cluster motion such as radial 
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pulsations (expansion and contraction), shear etc., that may have been originally 
present, have by now been completely damped out. The one mode that could have 
survived is rotation because angular momentum is conserved, although Weyman 
(1967) hnds that the rotation about three mutually perpendicular axes is consistent 
with zero to within 0.05 kms“^pc“^. 


(2): We assume a plausible shape for the isochrone of the Hyades main-sequence. However, 
this turns out to be not a very restrictive assumption since we £t for the amplitude 
and the slope of the main-sequence using the photometric data of the Hyades members 
itself. In principle, the shape of the isochrone could affect our selection of the 
Hyades members, although in practice, we hnd that the cluster membership remains 
unchanged for two different plausible isochrones. 


(3): We adopt the radial velocities collated from different sources in the literature as the 
astrometric radial velocities of the stars themselves, even though these might include 
contributions from non-astrometric sources such as convective and gravitational line 
shifts, atmospheric pulsations etc. (|Dravins, Larsson fc Nordlund 198^ ; |Nadeau,| 


19881) . However, most of the radial velocities used by us come from the work of 
Griffin et ah (1988), who calibrated the zero-point of their radial velocities using 
the radial velocities of asteroids and also applied a magnitude dependent correction 
to remove the contributions from convection in the stellar atmospheres (Gunn et ah 
1988, hereafter G88). In particular, 33 stars of our 43 Hyades members have their 
radial velocities measured by Griffin et ah (1988). 


(4): We assume that the parallaxes of the three Hyades spectroscopic binary systems 51 
Tauri, 70 Tauri and 6*^ Tauri derived by Torres, Stefanik & Latham (1997a, 1997b, 
1997c, hereafter T97a, T97b and T97c) are free from any systematic errors and that 
their quoted random errors are realistic. 


3. STATISTICAL PARALLAX AS A GENERALIZED MOVING CLUSTER 

METHOD 

In this section, we describe how the statistical parallax method is a generalized 
form of the classical moving cluster methods and provides a more precise distance to the 
Hyades cluster. First, we present the various equations describing the geometry of the 
cluster motion and show how the statistical parallax method subsumes the moving cluster 
methods. We then use some simple assumptions about the Hyades cluster to present a 
quantitative estimate of how much extra information about the distance is present in the 
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statistical parallax formalism compared to the moving cluster methods. However, while 
we present here a full statistical parallax solution (including the distance scale), the only 
real use that we make of this solution in the present paper is in selecting the cluster 
members. All the distance scale information in this solution is “over-ridden” by Hipparcos 
trigonometric parallaxes before comparing with the binary orbital parallaxes so that the 
latter can directly test the former. To make proper use of the statistical parallax Hyades 
distance measurement would require a careful examination of systematic errors and we 
defer this investigation to another paper. We note that once we select the members (and 
adopt Hipparcos parallaxes), the cluster space velocity does not depend on the details of 
the statistical parallax method at all. 

Consider a cluster at a distance d whose bulk velocity is V. If the radial velocity at an 
appropriately defined cluster center is K and the transverse velocity of the cluster in the 
plane of the sky is Vy, we have, 

Vr 

T’ 


= 


( 1 ) 


Vt = V - Kf, 


( 2 ) 


where is the proper motion vector of the cluster center in the plane of the sky. The 
difference between the transverse-velocity and the proper-motion vectors {6\t and 6fi) of 
the cluster center and those of the individual cluster members are then given by 


5Vt = -VrO, 


(3) 


5/x = 


'6Y2 

~ir 



(4) 


where 0 is the angular separation vector between the cluster center and the cluster 
member star in the plane of the sky, and we have assumed that |0| 1 (the small angle 

approximation), and {Sd/d) -C 1. This vector can be split into two components p|| and 
along the directions that are parallel and perpendicular respectively, to the proper motion 
vector (/x) of the cluster in the plane of the sky. Equation (^) can then be written in terms 
of these components as 

= - (^) 9i. (5) 

since = 0, by definition. Further, we also have 


SVr = (0-/x)d = 0||/i||d = 0||I/r, 


( 7 ) 
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where bVr is the difference between the radial velocities of the cluster member star and the 
cluster center, and Vt = iVrl. 


It is clear from the above equations that there are three independent measures of 
the distance to the cluster from equations (H), (j|), and (^. In the classical moving 
cluster method, the proper motions of the individual cluster members are used to derive a 
convergent point on the sky. This information is combined with an average radial velocity 
of the cluster center to derive its distance using equation (|]) ( [Boss 1908| ; |Hanson 1975| ; 
^chwan 199 1|) . Alternatively, if there are reliable radial velocities of the cluster members, 
equation (|^ can be used to estimate the cluster distance by making use of the average 
proper motion of the cluster center ( petweiler et al. 1984| ; G88). So far, equation has 
been used only in a restricted sense, in which the first term is used to derive a distance 
estimate, while the second term is simply neglected assuming that the spread in the 
distance of individual stars (5(i) adds to the uncertainty in this estimated distance (|Upton 
19701) . However, in this paper, we will use the photometric distance modulus to each star to 


estimate the quantity { 6 d/d), and hence reduce the uncertainty arising from the non-zero 
depth of the Hyades cluster. 


All three independent estimates of the cluster distance are naturally combined in 
the statistical parallax method. The resultant distance is then the weighted average of 
the individual distances from the three equations. Since these distance estimates are 
independent of each other, their variances add harmonically. The weight from each of these 
estimates is given by Wi = N{di/ai)^ where di and ai, {i = 1,2,3) are the distances and 
the errors in the distances from each of the three equations and N is the total number of 
members that are used to estimate the cluster distance. These weights are approximately 
given by 


ITi = N' 


(e^Vr, 


W 2 =N 


(dcr^)2 -h 0-2 
{OwVrf 


(d(T^)2 + a‘^ + {adlif 


W 3 = N 






( 8 ) 

(9) 

( 10 ) 


where and are the errors in the radial velocities and the proper motion respectively, cr^ 
is the uncertainty in the relative distance to individual cluster members, and a is the velocity 
dispersion of the cluster. The weight Wi corresponds to the classical convergent-point 
moving cluster method using individual proper motions [eq. (^], while IT 2 corresponds to 
the extension of this method using photometry to estimate the relative distances between 
the cluster members [eq. (H)]. The weight W 3 corresponds to the radial-velocity gradient 
method described by equation (|^ 
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For the purpose of illustration, we assume that for the Hyades cluster, a = 0.3 
kms“^, da^ = 0.3 kms“^, cr^ = 0.2 kms“^, = 0.6 kms“^, = (6*^^) = (6*^), and 
Vr = {5/3)Vt = 40 kms“^. This leads to hFi : hF 2 : = 1 : 0.33 : 0.5, showing that 

there is signihcant information about the distance in the two terms of equation (^. Here, 
we have assumed that ad = 1.5 pc, arising from an error of 0.06 mag in the photometric 
distance modulus of individual stars. Hence, the distance estimate using the statistical 
parallax method is more precise and accurate compared to that derived using the classical 
moving cluster methods alone. We note that in the absence of observational errors, the 
fractional accuracy in the cluster distance from the statistical parallax method using N 
cluster member stars is given by 


Ad 


= a 


( 21/2 + 1 /, 




- 1/2 


( 11 ) 


4. MEMBERSHIP AND COMMON CLUSTER MOTION 

We now hnd the bulk velocity of the Hyades cluster using the Hipparcos proper motions 
of the cluster members. However, a non-trivial problem here is the identihcation of the 
stars belonging to the cluster itself. We present our criteria for selecting the Hyades cluster 
members in §4.1 and describe our method of deriving its space velocity in §4.2. In §4.3, 
we compare our estimates of the cluster space velocity and distance modulus with previous 
determinations of these quantities in the literature. 


4.1. Cluster Membership 

We select a preliminary set of 75 Hyades cluster candidates from the list of 282 Hyades 
candidates in Table 2 of Perryman et ah (1998, hereafter P98). Our selection criteria for 
choosing these candidates are listed below in the order in which they are enforced. 

(1) : We reject all the stars that are either known or suspected to be binaries from earlier 

work. These are the stars with any alphabetical entry in at least one of the columns 
s, t or u of Table 2 of P98. This criterion rejects a total of 130 stars. 

(2) : All the candidate stars should have radial velocity measurements. We reject 14 stars 

under this criterion. We correct the raw radial velocity measurements of those stars 
measured by Griffin et al. (1988) (stars with an entry 1 in column r of Table 2 of P98) 
using the procedure described by equation (12) of G88, but accounting for a sign error 




([P98| ; R.Griffin 1998, private communication). Of the 43 cluster members used in this 
paper, we hnd that 33 stars have radial velocities measured by Griffin et ah (1988). 

(3) : We reject any star that is flagged as a variable in the Hipparcos catalog. This condition 

eliminates one star, with Hipparcos ID HIP 17962. 

(4) : We select only those stars that have ground based photometric measurements of both 

Vj and {B — V)j. We use the mean values of these quantities for the candidate stars 
from the GGPD photometric database of Mermilliod, Mermilliod & Hauck (1997), or 
from the Hipparcos catalog itself if the former are not available. 

(5) : The candidate stars should have colors in the range 0.4 < {B — V)j < 1.0. 

The last two photometric criteria (4) and (5) together eliminate another 62 stars, leaving us 
with 75 single star Hyades candidates in the color range 0.4 < [B — V)j < 1.0, all of which 
have ground based photometric data, radial velocity data and Hipparcos astrometric data. 

We derive the photometric distance modulus (m — M) to each of the 75 candidate stars 
in the color range 0.4 < {B — V)j < 1.0 by hnding the difference between the apparent 
magnitude of the star and the absolute magnitude for its color predicted by the isochrones 
of the Hyades main sequence. For the adopted range in {B — V)j color, the isochrones are 
reliable indicators of the distance modulus up to a possible global offset. Since the Hyades 
isochrones have not previously been determined to high precision, we apply our selection 
criteria using two distinct isochrones which, as we show below, span the range of the true 
isochrone. First, we use the isochrones adopted by PSSKH98 and we refer the reader to that 
paper for further details about the construction of the isochrones. We assume a metallicity 
of [Fe/H] = +0.14 and an age of 625 Myr for the Hyades (P98|) . We use the Yale color 
calibration ([Green 1988|) to transform the isochrones from the luminosity-temperature plane 
to the color-magnitude plane. Second, we use the color calibration proposed by Alonso, 
Arribas, & Mart Inez-Roger (1996) which predicts a different shape for the isochrone. These 
two isochrones have different zero points and color dependence with the result that if 
the isochrones are forced to coincide at {B — V)j = 0.4, they differ by about 0.3 mag at 
{B — V)j = 1.0. We assume that the true isochrone is in the general range of these two 
hducial isochrones and parametrize it by the function 

My{B -V) = Mv,Y.ie{B -V)+ A(m - M) + a[{B - V)j - 0.7] (12) 

where A(m — M) and a are parameters to be determined. These allow for both an offset in 
the zero point and a different slope for the color-magnitude relation. 

For each pair of values of A(m — M) and a, we determine the space velocities (Vj) of 
all these candidates using their photometric distance modulus, their radial velocities and 
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their proper motions from the Hipparcos catalog. We derive a best-ht mean velocity (V) 
of their centroid and reject the stars that are gross outliers from this mean cluster motion. 
We compute the quantity dehned as, 

N rp 

x" = E (Vi - v) C-‘ (v. - v) , ( 13 ) 

i=l 

where the summation is over all the N stars that remain at each iteration. The covariance 
matrix for star i includes contributions from the error in the photometric distance 
modulus, from the error in the radial velocity, from the error in the proper motion, and 
from the velocity dispersion of the cluster. The cluster velocity dispersion is between 0.2 
to 0.4 kms“^ (|G88| ; P98; D97) for plausible values of the cluster mass of about BOOM© to 
450Mq and half-mass radius of the cluster of about 4 to 5 pc (Pels, Oort fc Pels-Kluyver 

0.4 kms“^. We estimate the error in the distance modulus as arising solely due to the error 
in the {B — V)j color and assume an average slope of 6 for the isochrone to translate this 
to an error in [m — M). 

We iterate the procedure described above until there are no strong outliers and 
the velocities of all the remaining stars are consistent with a common cluster motion. 

At each iteration, we include only those stars with photometric distance modulus 
(m — M) < 4.5 -|- A(m — M), thus excluding stars that are located more than 30 pc beyond 
the centroid of the cluster. This procedure eliminates 13 stars from our original list of 75 
candidates, of which 12 are classihed as non-members by P98. The remaining one star 
(with Hipparcos ID HIP 20205) is a giant (spectral class G8HI according to the Hipparcos 
catalog) and we do not include it in our distance determination as its distance modulus 
cannot be estimated using an isochrone of the Hyades main sequence. We reject any star 
whose individual contribution to is greater than [Sx^/N) as an outlier. We derive our 
best estimates of a and A{m — M) from the remaining stars and use these values to estimate 
the cluster distance, its bulk velocity, and its velocity dispersion (a) using the statistical 
parallax algorithm. These two steps are described in detail in §4.2. The statistical parallax 
method also yields the extra error in the photometric distance modulus (e) which should be 
added in quadrature to the quoted errors so that the velocity dispersion tensor is isotropic. 
We repeat the cluster membership using these new values of a, A(m — M), a and e. We 
iterate the entire procedure until the cluster membership converges. 

We present the results of our iterative procedure in Table |l] in which we list the number 
of stars used at each iteration (Wtar), the total x^ cit the end of each iteration, the threshold 
value of the individual x^ of every star below which we accept the candidates as members 
(= 3x^/Astar), the number of stars that are rejected as non-members at each iteration 


1975; G88|, P98). In the hrst iteration of our membership selection, we assume a value of 
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(^reject,) and the number of rejected stars that are classihed as members by P98 (iVpgg). Of 
the remaining 62 Hyades candidates, 11 stars are classihed as non-members by both P98 
and by our algorithm. There are an additional 8 stars (with Hipparcos IDS HIP 20187, 
21788, 16908, 21112, 21267, 20491, 21066, 13806) that are classihed as members by P98, 
but which we reject as non-members because their individual contributions to the total 
are 235.2, 230.5,19.4,18.0,16.4,14.9,14.9, and 11.9 respectively, while the threshold value 
of {3x^/N) = 9.9 at the end of the last iteration. We note that even if we increase the 
threshold to (dy^/iV), only one additional star (HIP 13806) will be classihed as a member. 
On the other hand, the 5 highest individual contributions to of our member stars are 
9.3, 9.1, 7.0, 6.7 and 6.7 (corresponding to the stars HIP 19504, 20130, 22566, 20082 and 
20349). 

We select a total of 43 stars as Hyades cluster members at the end of the last iteration, 
all of which lie in a tight cluster around the mean cluster motion in velocity space. The best 
ht values of the parameters at the last iteration are a = 0.28 ±0.07, A{m — M) = 0.10 ±0.03 
mag, a = 320 ±39 ms“^, and e = 0.042 mag. Our membership selection procedure is robust 
to any changes in the absolute calibration of the isochrones since the relative distances 
between the cluster candidates are unahected by this. However, it is sensitive to the shape 
of the isochrones, although in practice, we hnd that the cluster membership is the same 
for the two hducial isochrones, despite the fact that they span a much larger range than is 
allowed by our hts for A(m — M) and a (see §4.2). 


4.2. Distance and Space velocity of Hyades 

We determine the common space velocity of the cluster from the velocities of all the 
cluster members selected by the procedure described in §4.1. We evaluate the y^ [as dehned 
in eq. (|^] in a dense grid of points in the space of the hve parameters namely, A(m — M), 
a and the three components of the bulk velocity (in equatorial coordinates) of the cluster 
(V). We ht this to a quadratic in the 5 parameters to hnd the best ht values and the 
covariance matrix of the parameters at the minimum of the y^ surface. We adopt this 
best-ht value of the bulk velocity as our initial guess of the cluster space velocity in the 
statistical parallax method. We estimate the components of the cluster space velocity in a 
coordinate system that is oriented such that one of the axes is along the radial direction of 
the centroid (W), another axis is along the direction perpendicular to the proper motion of 
the cluster in the plane of the sky (W) and the third axis is parallel to the proper motion 
of the cluster in the plane of the sky (V||). By dehnition, W is the radial velocity of the 
cluster, V\\ is its velocity in the plane of the sky, and V± is zero. 
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To compute the photometric distance to each star, we fix the slope-correction (a) 
and the zero point offset [A(m — M)] at their best-fit values derived as described above. 

In the modern version of the statistical parallax method as described by Popowski & 
Gould (1998), one uses the maximum likelihood procedure to simultaneously solve for ten 
different parameters viz., the distance scaling factor relative to a fiducial distance scale 
( 77 ), the three components of the bnlk velocity of the cluster (G-, V± and V)[) and the six 
independent components of the second moments of its velocity distribution - the three 
diagonal terms corresponding to the sqnare of the velocity dispersion in the three directions 
(a^, a]_ and ajj) and the three nnique off-diagonal terms (cr^j_, cr^w and We assnme an 

isotropic velocity dispersion tensor of the Hyades with the result that the three independent 
off-diagonal terms are constrained to be zero. For an assumed level of errors in the radial 
velocities, the proper motions, and the distance to individual stars, the statistical parallax 
method derives a maximum likelihood estimate of the cluster velocity dispersions in the 
three mutnally perpendicular directions. However, the errors in the distance to each star 
affect only the parallel dispersion a\\ while the estimates of the velocity dispersions in the 
radial and the perpendicular directions (uj. and a± respectively) are independent of the 
distance errors. Therefore, we begin by constraining the velocity dispersion in these two 
directions to have the same valne. The velocity dispersion in the parallel direction (an) 
now includes contributions from both the intrinsic velocity dispersion of the cluster and the 
dispersion arising from a possibly wrong estimate of the distance errors. Hence, we add an 
extra error (e) in qnadratnre to onr original errors (as listed in the photometric sources) in 
the photometric distance modulus so that the velocity dispersion in the parallel direction 
(a||) becomes equal to the velocity dispersions in the other two directions. We adopt this 
velocity dispersion as the velocity dispersion of the cluster. 

At the end of the last iteration, the sample contains 43 Hyades cluster members, the 
velocity dispersion is eqnal to 320 ±39 ms“^, and we need to add an extra error of e = 0.042 
mag (in quadrature) in the photometric distance modulus to enforce an isotropic velocity 
dispersion tensor. We adopt this value of the velocity dispersion and e in the remainder 
of this paper. The statistical parallax method finds the maximum likelihood solution of 
the five independent parameters namely, p = (77, W, V±, V\\, a^) subject to the constraint of 
an isotropic velocity dispersion tensor. This solution, using only the relative photometric 
distances to individnal Hyades cluster stars is given by p(stat) = [1.0027 ± 0.0143, 39.50 ± 
0.06 kms-^ 0.00 ± 0.07 kms-^ 23.82 ±0.34 kms-^ 0.1034 ± 0.0251 (kms-i)^], where the 
velocity dispersion is, as expected, a^ = 320 ± 39 ms“^. 

The statistical parallax measurement of the distance scale (77phot) is of interest in its 
own right. However, a proper interpretation of this measurement reqnires a thorongh 
investigation of the systematic errors which is beyond the scope of this paper. Here, our 
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primary goal is to test the Hipparcos trigonometric parallax distance scale toward the 
Hyades. To this end, we will therefore renormalize rj to force agreement (on average) 
between the distances to cluster members as determined in the renormalized solution and 
from Hipparcos parallaxes. This will have the effect of rescaling V\\ by the same factor, but 
will have no effect on our estimate of 14 or Vj_. The resulting space velocity of the cluster 
( 14 , V±, 1 ^ 1 ) can then be used to measure the distances to individual Hyades binary systems 
from their proper motions and thus test the Hipparcos trigonometric-parallax Hyades 
distance scale against the orbital parallaxes of these binaries. 


We begin by computing for each cluster member, the quantity 


Vi = 


^photd 


where 


^7r,Hip,2 


(1 -h , 


Xi = 


(14) 


(15) 


^phot,2 

Since the errors in the Hipparcos parallaxes are approximately Gaussian distributed 
( lArenou et al. 19^ ), the quadratic correction term x^ is required to ensure that the two 
sides of equation (^) have the same mean value (|Lutz fc Kelker 1975| ; [Smith fc Eichhorii 
1996 ). The error in rji is given by 


a 


r],i 



(16) 


We hnd that the mean value of p for all the cluster members is given by r/Hip = 1.0178±0.0081. 
This is consistent with the value of Pphot = 1.0027 ± 0.0181 derived using the Hipparcos 
proper motions alone. The difference between these two distance scaling factors is given by 


Vphot - Vmp = -0.0151 ± 0.0164, (17) 

and hence is consistent with zero, thus providing a semi-independent check of the 
self-consistency of Hipparcos astrometry. Although P98 reached a similar conclusion, 
our result is based on improved parallax estimates, as the statistical parallax solution 
also includes the photometric distance information to individual stars. This consistency 
between the scaling factors would permit us, if we desired, to combine the two 
determinations to hnd a single best estimate of p. However, our aim here is to test 
the Hipparcos parallaxes alone. Therefore, we simply constrain the overall solution to 
have 7 ] = pHip = 1.0178 ± 0.0081. We then hnd that the space velocity of the cluster is 
(14, hj_, V||) = (39.51 ± 0.06 kms“^, 0.00 ± 0.07 kms“^, 24.17 ± 0.22 kms“^) and the matrix 
of correlation coefficients is 

/ 1.0000 -0.0045 0.0161 \ 

-0.0045 1.0000 -0.1392 

^ 0.0161 -0.1392 1.0000 y 


( 18 ) 
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The two primary effects of using the Hipparcos parallaxes to hx the distance scale to the 
Hyades members are: 

(1) : It increases the value of V\\ by 1.5% i.e., almost exactly equal to the increase in the 

value of rj itself. 

(2) : It signihcantly reduces the error in V\\. 

However, neither nor V±, nor their errors are signihcantly affected. The fact that 
the maximum likelihood value of rj is greater than one means that the Hipparcos 
parallax distance scale to the Hyades is larger than our hducial scale by a factor 
(?7 — 1) = 1.78%. The space velocity of the cluster in equatorial coordinates is 
{Vx,Vy,Vz) = (—5.70 ± 0.20,45.62 ± 0.11, 5.65 ± 0.08) kms“^ and the matrix of correlation 
coefficients is 

/ 1.0000 -0.8134 0.5042 \ 

-0.8134 1.0000 -0.4964 . (19) 

^ 0.5042 -0.4964 1.0000 y 

We will use this estimate of the bulk velocity of the cluster in the remainder of the paper to 
predict the parallaxes from the Hipparcos proper motions of individual stars. In Galactic 
coordinates, this velocity is given by {U, V, W) = (—42.27±0.08, —18.89±0.18, —1.51T0.14) 
km s“^. 

We show the velocities of the Hyades cluster candidates in Figure 1. The hrst three 
panels (a)-(c) show the velocities computed using the photometric distance modulus 
to each star (normalized to the Hipparcos trigonometric parallax distance scale) and 
the Hipparcos proper motions, while panel (d) shows the velocities computed using the 
Hipparcos parallaxes and proper motions. The crosses show the velocity components of 
the 43 cluster members, while the open circles represent the velocity components of the 8 
stars that are classihed as members by P98, but rejected as non-members by our algorithm. 
The smaller scatter in velocities of the members in panel (c) compared to that in panel (d) 
shows that the photometric distance moduli lead to a much tighter core in velocity space, 
and hence a cleaner separation between the members and the non-members compared to 
using distances inferred from Hipparcos parallaxes. The solid circle in all the panels shows 
the bulk velocity of the cluster. For this space velocity of the cluster, the total is 143 for 
43 stars (corresponding to 126 degrees of freedom) demonstrating that our estimates of the 
errors for the various quantities and of the cluster velocity dispersion are reasonable. The 
centroid is at a distance of |r| = 46.61 ± 0.38 pc (corresponding to a distance modulus of 
3.34 ± 0.02), and its equatorial coordinates are a = 04^26'"32^, <5 = 17°13.'3 (2000). This is 



also the direction of the radial velocity of the cluster center, i.e, the direction of 14-. The 
motion of the cluster in the plane of the sky is towards the direction 105°11.'0 East of North. 


4.3. Comparison with previous estimates 

We compare our estimates of the distance modulus to Hyades centroid, its bulk space 
velocity, and its velocity dispersion, with previous determinations of these quantities in the 
literature. 

Our distance modulus to Hyades of (m — M) = 3.34 ± 0.02 mag agrees very well 
with the value of 3.33 ± 0.01 mag obtained by P98 using the stars located within 10 pc of 
the cluster center, and the value of 3.34 ± 0.04 mag determined by PSSKH98 using the 
main-sequence htting technique. We would like to emphasize here that while the P98 value 
refers to the distance modulus of the center of mass of the Hyades, our estimate refers to 
the centroid of the 43 Hyades members selected in §4.1. However, the estimate of the bulk 
velocity does not depend on the choice of the subsample of the cluster and our results can 
therefore be directly compared to those of P98 and D97. 

The space velocity of the Hyades cluster has been recently determined by 
P98 and by D97 using the Hipparcos astrometric data of the Hyades members. 

We compare these velocity estimates with our estimate of the bulk velocity of 
the cluster in the coordinate system centered on the centroid of Hyades, viz., 

(W, 14, 4|) = (39.51 ± 0.06 kms-^ 0.00 ± 0.07 kms-^ 24.17 ± 0.22 hms"^). In this 
coordinate frame, the cluster velocity determined by P98 is (38.82,0.03,24.55) kms“^, 
while that determined by D97 is (39.60,0.00,24.65) kms“^. The errors in these velocity 
components are unlikely to exceed 0.2 kms“^. Our estimate of the radial velocity of the 
Hyades centroid is signihcantly larger than that of P98, but agrees very well with the value 
of D97. However, the first two panels in Figure 1 show that the radial velocities of 5 of the 8 
stars that are classihed as members by P98, but rejected as non-members by our algorithm 
are signihcantly (more than 1 kms“^) smaller than the mean radial velocity of the cluster, 
and including them as members will systematically reduce the mean cluster radial velocity. 
These 5 stars are clearly well outside the tight cluster in the radial velocity component 
(W), have high values of and are therefore unlikely to be Hyades members given that 
the cluster velocity dispersion is only 320 ± 39 ms“^. Although the radial velocities of 
the other 3 stars (HIP 13806, 21066 and 21112) are consistent with them being Hyades 
members, they have systematically high values of the parallel velocity component (Vji) and 
are therefore classihed as non-members. However, the star HIP 20350 whose radial velocity 
is about 41.44 kms“^ is accepted as a Hyades member because of the large error in its radial 
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-4 -2 0 2 4 -4 -2 0 2 4 

Vj^ [kms“*] Vj^ [kms“*] 


Fig. 1.— Velocities of the Hyades cluster candidates from the Hipparcos proper motions. 
The velocities in panels (a)-(c) are computed using the photometric distance to each star 
(normalized to the Hipparcos trigonometric parallax distance scale), while those in panel 
(d) are computed using the distance inferred directly from the Hipparcos parallax, (a) V- 
and V± (b) V- and V\\ (c) V± and V\\ and (d) V± and V\\. The component K is the velocity 
component in the radial direction of the centroid, while V± and V\\ represent the velocity 
components perpendicular and parallel to the direction of the proper motion of the cluster 
in the plane of the sky. The crosses are the velocity components of the 43 cluster members 
while the open circles are those of the 8 stars that are classihed as members by P98, but 
rejected as non-members by our algorithm. The solid circle in each panel shows the mean 
motion of the cluster. 
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velocity (cr^ = 2.4 kms ^). Hence, the difference between our estimate of the radial velocity 
and that of P98 can be attributed to the differences in the selection of Hyades members. 


All the estimates of the velocity component in the perpendicular direction (Vj_) are 
consistent with each other, and are individually consistent with the expected value of 
zero. However, our estimate of the parallel velocity component (V|| = 24.17 ± 0.22) is 
inconsistent at the 1.7cr level with the value derived by P98, and at the 2.2cr level with 
that derived by D97. The cluster velocity dispersion of 320 ± 39 ms“^ computed using our 
statistical parallax algorithm agrees well with the value of 300 ms“^ determined by P98 
and acceptably with the value of 250 ± 40 ms“^ determined by D97, although our value is 
computed excluding all the binary systems that dominate the mass distribution near the 
central regions of the cluster (P98, [Pels, Uort fc Pels-Kluyver 1975|; [lerlevich 1987|). 


5. PARALLAX FROM PROPER MOTION 


We adopt the cluster space velocity derived in the previous section to predict the 
parallaxes of the 43 member stars using their proper motions from the Hipparcos catalog. 
The parallax of any cluster member that has the same space velocity as the cluster is given 
by 


^pm,2 


((v.)ilc.-V,) 

((V.).|C-‘|(V.)i) 


( 20 ) 


where (Vt), = Vc - (fi ■ Vc) f j is the transverse velocity of the cluster in the plane of the sky 
at the position of the star i, is its proper motion from the Hipparcos catalog and C,, the 
covariance matrix, is the sum of the proper motion error tensor of star i and the velocity 
dispersion tensor divided by the square of the distance. We have employed Dirac notation. 


{X|0|Z) = yjfjC>.,Z,, (21) 

i,j 

1 /2 

The error in TTp^,* is equal to ((Vt)i|C"^|(Vt)i^ . 

We show the difference between the photometric parallax (TTphot) and the parallax 
determined assuming a common space velocity for all the cluster members (vTpjn) in Figure 2. 
We have scaled all the photometric distances by the factor r^Hip = 1-0178 so that the mean 
value of this difference should be equal to zero. The horizontal error bars show the error in 
photometric parallax, while the vertical error bars show the uncertainty in TTp^ alone. Thus, 
the total error in the quantity (TTphot — T^pm) is the sum in quadrature of the two error bars 
displayed in Figure 2. We find that the mean offset between the two parallaxes is given by 


(TTphot — TTpm) = —0.009 ± 0.099 mas. 


(22) 
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while the of this difference is 29.2 for a total of 43 stars. This shows that there are no 
internal inconsistencies in our method of predicting the parallaxes from the proper motions. 


6. PARALLAXES OF BINARY SYSTEMS 

We predict the parallaxes of the 3 binary systems 51 Tauri (HIP 20087), 70 Tauri 
(HIP 20661) and 6 ^ Tauri (HIP 20894) using the cluster space velocity (Vc) determined 
above and their individual Hipparcos proper motions. The full orbital solutions of these 3 
spectroscopic binary systems have been derived by T97a, T97b and T97c respectively. The 
Hipparcos proper motions refer to the center of mass for HIP 20087, while they refer to the 
motion of the photocenter for the other two systems. HIP 20894 whose semi-major axis is 
less than 0."1 is listed as a variable single star in the Hipparcos catalog. For HIP 20661 and 
HIP 20894, we compute the difference between the proper motion of the center of mass and 
the proper motion of the photocenter using the spectroscopic-astrometric orbital solutions 
for these 2 binary systems provided by T97b and T97c respectively. In Table 2, we list the 
proper motions of the center of mass of all the 3 systems, their parallaxes from Hipparcos, 
their parallaxes from their proper motions, and their orbital parallaxes. There are two 
important features in the errors of the different parallaxes in Table ^ 

(1) : The error in the parallaxes determined from the individual proper motions of the 

binary systems and the common space velocity of the cluster [(J7r(pm) in column 7 of 

Table 2] is almost a factor of three smaller than the error in the Hipparcos parallaxes. 

(2) : The errors in the orbital parallaxes of these binary systems [(j^(orb) in column 9 of 

Table 2] are about twice as large as the errors in the proper-motion parallaxes. 

The mean difference (weighted by the inverse square errors) between the proper-motion 
parallaxes and the orbital parallaxes for these three binary systems is, 

(TTpm — TTorb) = 0.52 ± 0.43 mas. (23) 

The error in this difference is dominated by the error in the orbital parallaxes. From Table 
1 , we find that the error in the individual proper-motion parallaxes are each of the order of 
0.3 to 0.35 mas. Hence, it appears possible in principle to detect any systematic errors in 
the Hipparcos parallaxes at the level of agys = [X] l/o'^(pKi)] = 0.2 mas if the orbital 

parallax errors could be reduced below this level. However, there is another source of error 
in determining the proper-motion parallaxes. This arises from the error in the component of 


Table 1: Cluster membership iterations 


Iteration 

star 


(3xVA^star) 

^reject 

Npgs 

1 

62 

23282 

1127 

6 

0 

2 

56 

3947 

212 

5 

0 

3 

51 

682.3 

40.14 

2 

2 

4 

49 

235.2 

14.4 

2 

2 

5 

46 

182.2 

11.9 

2 

2 

6 

44 

153.8 

10.5 

1 

1 

7 

43 

142.5 

9.9 

0 

0 


Table 2: Astrometry of the 3 spectroscopic binary systems with orbital parallaxes. 


HIP ID 

/iQ,cos((5) 
(mas yr“^) 

(mas yr“^) 

TJ'ffip 

(mas) 

a^(Hip) 

(mas) 

^pm 

(mas) 

a^(pm) 

(mas) 

^orb 

(mas) 

(T7r(orb) 

(mas) 

20087 

96.42 

-33.92 

18.25 

0.82 

18.45 

0.30 

17.92 

0.58 

20661 

104.97 

-26.67 

21.47 

0.97 

21.35 

0.36 

21.44 

0.67 

20894 

108.80 

-26.35 

21.89 

0.83 

22.49 

0.37 

21.22 

0.76 
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Fig. 2.— Difference between the photometric parallax (TTphot) and the parallax predicted 
assuming a common space velocity for the cluster members (TTpm). All the photometric 
distances have been scaled by the quantity T^Hip = 1-0178 so that the mean value of this 
difference should be zero. The horizontal error bars show the error in TTphot, while the 
vertical error bars show the uncertainty in TTpm alone, i.e, the error in (TTphot — TTpm) is the 
quadrature sum of the two error bars. 
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the cluster space velocity itself in the direction of its proper motion in the plane of the sky, 
i.e, in the component V^|. This error is equal to {nuya} (<^V||/^|) = 0-19 i^^s where we have 
again assumed that the mean parallax of the Hyades cluster is 21.5 mas. Adding this error 
in quadrature to the errors determined above, we hnd that (TTpm — TTorb) = 0.52 ± 0.47 mas. 
The irreducible error in this method Usys is now 0.28 mas, still considerably smaller than 
the errors in the currently available orbital parallaxes. We note that this analysis assumes 
that there are no systematic errors in the orbital parallaxes of the binary systems. 

P98 also concluded that the Hipparcos parallaxes of these binary systems and their 
parallaxes from their orbital solutions are consistent. However, the accuracy of their 
comparison is limited by the errors in the Hipparcos parallaxes. Thus, we hnd that the 
mean difference (weighted by the inverse square errors) between the Hipparcos parallaxes 
and the orbital parallaxes for these three binary systems is 

('TJ'Hip — TJ'orb) = 0.35 ± 0.63 mas. (24) 

The error in this difference is dominated by the error in the Hipparcos parallaxes and is 
therefore irreducible in the future. This shows that the technique proposed in this paper 
allows one to check the level of systematic errors in the Hipparcos parallaxes at more than 
twice the precision than is possible using the straightforward comparison of Hipparcos and 
orbital parallaxes. 


7. CONCLUSIONS 

Our main conclusions are as follows: 

(1) : When the distance scale to the Hyades cluster is hxed by Hipparcos parallaxes, the bulk 

velocity in equatorial coordinates is (14, Vy, U) = (—5.70±0.20,45.62±0.11, 5.65±0.08) 
kms“^, its velocity dispersion is 320 ± 39 ms“^ and the distance modulus of the 
centroid of our sample of 43 members is {m — M) = 3.34 ± 0.02. This distance 
modulus agrees with that determined by both P98 and PSSKH98 and the velocity 
dispersion is consistent with the estimates of both P98 and D97. However, the 
difference between our estimate of the bulk velocity and that of P98 arises from our 
more stringent selection of Hyades members. 

(2) : The Hipparcos parallaxes of the three Hyades binary systems are consistent at the la 

level with the parallaxes from their orbital solutions. Hence, the systematic error in 
the Hipparcos parallaxes towards the Hyades cluster is less than 0.47 mas. This is 
already more precise (by a factor of almost 1.4) than a straightforward comparison 
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between the Hipparcos and orbital parallaxes performed by P98, in which the precision 
is limited by the Hipparcos parallax errors, and hence is irreducible in the future. 

(3): The test proposed in this paper can, in principle, detect any systematic error greater 
than 0.3 mas in the Hipparcos parallaxes towards the Hyades cluster. The dominant 
factor that currently limits a check at this level is the “large” errors in the orbital 
parallaxes of the binary systems. 

It follows from the last two points that a more accurate estimate of the binary orbital 
parallaxes would enable a better determination of the systematic errors (or the lack thereof) 
in the Hipparcos astrometry towards the Hyades cluster. 


This work was supported in part by the grant AST 97-27520 from the NSF. We thank 
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We are grateful to Bob Hanson for his perspicuous comments on an earlier draft of 
this paper. We also thank the referee for a detailed report that contained many useful 
suggestions. 


REFERENCES 

Arenou, F., Lindegren, L., Froeschle, M., Gomez, A.E., Turon, C., Perryman, M. A.C. & 
Wielen, R. 1995, A&A, 304, 52 

Arenou, F., Mignard, F. & Palasi, J. 1997, Chapter 20, Verihcation of Parallaxes, in the 
Hipparcos and Tycho Catalogues, Vol. 3, Construction of the Hipparcos Catalogue, 
ESA SP-1200 

Alonso, A., Arribas, S., & Mart Inez-Roger, C. 1996, A&A313, 873 
Boss, L., 1908, AJ, 26, 31 

Dravins, D., Larsson, B., & Nordlund, A., 1986, A&A, 158, 83 

Dravins, D., Lindegren, L., Madsen, S., Holmberg, J. 1997, in Hipparcos Venice ’97, eds. B. 
Battrick & M. A. C. Perryman (ESA), 643 (D97) 

ESA 1997, The Hipparcos and Tycho Catalogues, ESA SP-1200 

Detweiler, H.L., Yoss, K.M., Radick, R.R. & Becker, S.A., 1984, AJ, 89, 1038 



- 22 - 

Hawley, S.L., Jeffreys, W.H., Barnes, T.G.III, & Wan, L., 1986, ApJ, 302, 626 
Hanson, R.B., 1975, AJ, 80, 5 

Green, E.M. 1988, in Galibration of Stellar Ages, ed. A.G.D. Phillip (Schenectady, 
NY:Davis), 81 

Griffin, R.F., Gunn, J.E., Zimmerman, B.A. & Griffin, R.E.M., 1988, AJ, 96, 172 
Gunn, J.E., Griffin, R.F., Griffin, R.E.M. & Zimmerman, B.A., 1988, AJ, 96, 198 (G88) 
Lutz, T.E. & Kelker, D.H., 1973, PASP, 85, 573 
Mermilliod J.-G., Mermilliod M., Hauck B. 1997, A&AS, 124, 349 

Mermilliod, J.-G., Turon, G., Robichon, N., Arenou, F.,& Lebreton, Y. 1997, in Hipparcos 
Venice ’97, eds. B. Battrick & M. A. G. Perryman (ESA), 643 

Nadeau, D., 1988, ApJ, 325, 480 

Popowski, P. & Gould, A., 1998, ApJ, 506, 000 

Pels, G., Oort, J.H., & Pels-Kluyver, H.A., 1975, A&A, 43, 423 

Perryman, M. A. G., Brown, A.G.A., Lebreton, Y., Gomez, A., Turon, G., Gayrel de 

Strobel, G., Mermilliod, J.-G., Robichon, N., Kovalevsky, J., & Grifo, F. 1998, A&A, 
331, 81 (P98) 

Pinsonneault, M., H., Stauffer, J., Soderblom, D.R., King, J.R. & Hanson, R.B., 1998, ApJ, 
in press (PSSKH98) 

Robichon, N., Arenou, F., Turon, G., Mermilliod, J.-G.,& Lebreton, Y. 1997, in Hipparcos 
Venice ’97, eds. B. Battrick & M. A. G. Perryman (ESA), 567 

Schwan, H., 1991, A&A, 243, 386 

Smith, H. & Eichhorn, H., 1996, MNRAS, 281, 211 

Strugnell, P., Reid, N., & Murray, G.A., 1986, MNRAS, 220, 413 

Terlevich, E., 1987, MNRAS, 224, 193 

Torres, G., Stefanik, R.P. & Latham, D.W., 1997, ApJ, 474, 256 (T97a) 

Torres, G., Stefanik, R.P. & Latham, D.W., 1997, ApJ, 479, 268 (T97b) 



Torres, G., Stefanik, R.P. & Latham, D.W., 1997, ApJ, 485, 167 (T97c) 
Upton, E.K.L., 1970, AJ, 75, 1097 
Weyman, P., 1967, PASP, 79, 156 


This preprint was prepared with the AAS MUX macros v4.0. 



